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Abstract: We present a design methodology and analysis 
of a cavity optomechanical system in which a localized GHz 
frequency mechanical mode of a nanobeam resonator is 
^ evanescently coupled to a high quality factor {Q > 10^) optical 

mode of a separate nanobeam optical cavity. Using separate 
nanobeams provides flexibility, enabling the independent 
^1 design and optimization of the optics and mechanics of the 

r-| ■ system. In addition, the small gap (^25 nm) between the 

^ O-t two resonators gives rise to a slot mode effect that enables 

a large zero-point optomechanical coupling strength to be 
achieved, with g/ln > 300 kHz in a Si3N4 system at 980 nm 
and g/ln ^ 900 kHz in a Si system at 1550 nm. The fact 
that large coupling strengths to GHz mechanical oscillators 
can be achieved in Si3N4 is important, as this material has a 
broad optical transparency window, which allows operation 
throughout the visible and near-infrared. As an application 
of this platform, we consider wide -band optical frequency 
conversion between 1300 nm and 980 nm, using two optical 
nanobeam cavities coupled on either side to the breathing mode 
of a mechanical nanobeam resonator. 
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1. Introduction 

Recent experiments using radiation pressure effects in cavity optomechani- 
cal systems [1], such as electromagnetically-induced transparency [2, 3] and 
ground state cooling [4, 5], as well as theoretical proposals for using such op- 
tomechanical systems in quantum information processing applications [6, 7, 8, 
9, 10], generally rely upon the achievement of large zero-point optomechanical 
coupling rates and low optical and mechanical dissipation rates. They also typ- 
ically require operation in the resolved sideband limit in which the mechanical 
resonance frequency must exceed the optical dissipation rate. One-dimensional 
optomechanical crystals [11, 3, 8], in which a periodic patterning is applied to 
a doubly clamped nanobeam in order to simultaneously localize GHz phonons 
and near-infrared (200 THz) photons, are one promising architecture for future 
work in quantum cavity optomechanics. 

While the use of a single nanobeam in such work has many advantages, sep- 
aration of the optical cavity and mechanical resonator can afford greater lev- 
els of flexibility. This can be useful in applications in which multiple optical 
modes are coupled to the same mechanical resonance, for example, to achieve 
wavelength conversion [6, 7, 12, 13] or in converting traveling wave phonons 
to traveling wave photons [6]. One disadvantage of separating the mechani- 
cal resonator and optical cavity is that the optomechanical coupling rate gen- 
erally tends to be smaller, given the reduced overlap between the mechanical 
and optical modes. To overcome this, one can focus on optical modes formed 
within the gap between the optical and mechanical resonator. Such air- slot type 
modes [14], which display a high electric field concentration in the air gap, 
have been used to achieve large optomechanical couplings [15, 16, 17, 18], and 
sub- wavelength effective optomechanical lengths [19]. In these demonstrations, 
however, the mechanical resonators and optical cavities were not separated, and 
mechanical resonance frequencies were typically < 150 MHz. For operation in 
the resolved sideband regime, higher frequencies are desirable, given that high- 
Q cavities in materials like Si and Si3N4 typically have optical dissipation rates 
that are > 100 MHz. 



Here, we present the design of a cavity optomechanical system in which a 
GHz frequency localized mode in a mechanical nanobeam resonator is coupled 
to a localized 300 THz optical mode in a separate optical nanobeam resonator 
(Fig. 1). We show that the slot mode effect can increase optomechanical cou- 
pling rates with respect to those found in a single nanobeam optomechanical 
system, and consider applications of this approach to instances in which the 
separated mechanical and optical resonators are particularly advantageous, such 
as wide-band optical wavelength conversion [6]. We consider both Si3N4 and Si 
systems, although our primary focus is on Si3N4. In particular, the Si3N4 system 
has many potential advantages for applications in quantum cavity optomechan- 
ics: optical 2 > 10^ has been demonstrated over a wide range of wavelengths in 
the visible and near-infrared in Si3N4 [20, 21, 22]. This can allow compatibility 
with quantum systems based on trapped atoms, ions, and semiconductor quan- 
tum dots, whose relevant optical transitions are often below 1000 nm (where Si 
is opaque). Furthermore, these stoichiometric Si3N4 films under tensile stress 
have been shown to support high mechanical quality factors (^ 10^), both in 
commercial thin films [23] and in nanofabricated devices [24, 25, 26]. Finally, 
as Si3N4 does not exhibit the two-photon absorption and subsequent free-carrier 
absorption and dispersion that silicon does [27], higher optical powers are us- 
able in this system, potentially increasing the range of pump-enhanced optome- 
chanical coupling values that can be achieved. Our geometry is predicted to 
support GHz mechanical frequencies, optical Q > 10^, and a zero-point op- 
tomechanical coupling strength g/2;r > 300 kHz. Such values should in princi- 
ple enable strong radiation pressure driven phenomena. 

The proposed geometry, shown in Fig. 1, consists of two suspended paral- 
lel dielectric beams of refractive index n^, widths Wq and Wm and thickness 
t, separated by a gap of width Wgap. Both beams are etched through to form 
quasi-periodic ID structures in the z-direction. The "optical" beam's ID lattice 
is modulated so that a quasi-TE (Ey = across the center of the Si3N4 layer at 
y = 0) confined optical mode is supported at an optical frequency /„. Optical 
confinement along the beam is provided by a ID photonic bandgap, while total 
internal reflection provides confinement in the x and y directions [11, 28, 29]. 
With a sufficiently small spacing Wgap, the quasi-TE mode consists of an air- 
gap resonance as in [14], with a strongly concentrated, almost completely x- 
oriented electric field between the two beams. The "mechanical" beam is de- 
signed to support a mechanical resonance of frequency /,„, confined in z, that 
displaces the boundary B facing the optical resonator, thereby closing the gap 
and thus shifting the optical resonance frequency fo. It is important to stress that 
the optical cavity is actually formed by the "optical" beam and the "mechani- 
cal" beam's closest arm (indicated in Fig. 1), so that optical resonances depend 
strongly on the gap width. Nevertheless, optical and mechanical resonances can 



Fig. 1 . Schematic of double beam optomechanical resonator. 

be tuned with great independence, as we show in the following sections. Sec- 
tions 2 and 3 present background information on the simulations performed 
and consider the case of a single Si3N4 nanobeam optomechanical cavity, as a 
baseline reference. Section 4 describes the design of the slot-mode-coupled op- 
tical and mechanical nanobeam resonators, emphasizing the design approach 
for achieving high-Q optical modes, GHz mechanical modes, and large op- 
tomechanical coupling. Finally, section 5 considers a few applications of this 
platform. 

2. Background 

The shift in the frequency fo of a particular optical resonance due to displace- 
ment of the nanostructure boundaries produced by a mechanical resonance at 
frequency /„ is quantified by the optomechanical coupling ^qm = dOo/dx = 
COo/Lqm', here, x is the cavity boundary displacement, cOo = 2nfo and Lqm is an 
effective optomechanical interaction length [30]. The effective length Lqm can 
be estimated via the perturbative expression [31] 



2jdV£\E\ 




/JA(|Q-n|)| 


fA£ E|| 


2 


-A(e-i) 


|D±I 


1^) 



(1) 



Here, E and D are the modal electric and electric displacement fields, respec- 
tively, Ae = £diei. - £air, A(e" ^ ) = E^^^j - E^-J: , and £diei.,air are the permittivities 
of the nanobeam material and air. The mass displacement due to the mechanical 
resonance is given by Q, and the normal surface displacement at the structure 
boundaries is |Q ■ n|, where n is the surface normal. The integral in the denom- 
inator is performed over the entire surface of the nano structure. 

The optomechanical coupling ^OM can be converted into a pure coupling 
rate g between the optical and mechanical resonances: g = x^pf ■ gou, where 
x^pf = ^Jh/lmcOm is the zero point fluctuation amplitude for mechanical dis- 
placement [6] and m is the motional mass of the mechanical resonance at fre- 
quency (Dm- The motional mass can be obtained from the displacement Q and 

the nanobeam material density p by m = p j dV \ ^}'^\Q\\ j ■ 



We point out that eq.(l) only quantifies frequency changes due to boundary 
distortions. An additional, potentially important contribution to the optome- 
chanical coupling comes from the photo-elastic effect, which corresponds to 
local, stress-induced changes of the refractive index. This contribution can be 
quite significant, and was in fact recently been exploited to demonstrate, in 
combination with the moving-boundary contribution, optomechanical coupling 
rates in excess of 1 MHz in silicon nanobeam optomechanical crystals [32]. 
Giant enhancement of Brillouin scattering in suspended silicon waveguides, 
resulting from both photo-elastic and radiation pressure optomechanical cou- 
pling, has also been predicted [33]. In the present work, only moving-boundary 
contributions were taken into account. 

Optical and mechanical resonator modes are obtained using the finite element 
method, as described in [15], and are used in the expressions above to yield the 
optomechanical coupling rates. 

3. Si3N4 nanobeam 

To gain some initial insight into the optomechanical rates achievable with a 
low index contrast nanobeam with 'standard' (i.e., single nanobeam) geometry, 
we start by investigating an optical cavity composed of a ID array of air holes 
etched on a Si3N4 nanobeam, as show in Fig. 2(a) [34]. We seek to achieve the 
largest possible optomechanical coupling rate g and optomechanical coupling, 
^OM- Our cavity was designed to support the optical mode at a wavelength 
A = 960 nm shown in Fig. 2(b). The 980 nm wavelength band is technologi- 
cally important due to the availability of triggered single photon sources based 
on InAs quantum dots [35] and silicon based single photon detectors which op- 
erate within this range [36]. Together with the optical mode, this cavity supports 
confined mechanical modes co-located with the optical resonance. The funda- 
mental breathing mechanical mode (FBM) at / = 3.2 GHz, shown in Fig. 2(e), 




with motional mass m = 570 fg, displays the highest optomechanical coupling 
rate g to the optical resonance in Fig. 2(b). 

To reach this design, the procedure outlined in Section 4.1 was initially 
employed to obtain a high Q optical mode. The resulting cavity displayed a 
quadratic variation of the lattice constant (the spacing between holes) from the 
center to the edges [28], and supported, together with the optical mode, a lo- 
calized breathing mechanical mode. The cavity geometry was then optimized 
via a nonlinear minimization code to yield a maximized optomechanical cou- 
pling rate g/27t = 133.6 kHz (Lqm = 5.1 /im), while maintaining a high opti- 
cal quality factor Q > 1 x 10^. The optimization was realized by allowing the 
aspect ratio of the elliptical holes (vertical axis over horizontal axis) to vary 
quadratically from the cavity center, then searching a 2D space of aspect ratios 
at the center and at the edge of the cavity. Reducing the aspect ratio at the cavity 
center proved to be beneficial in reducing the effective length Lqm, by reducing 
the mechanical mode volume and promoting better overlap between the opti- 
cal and mechanical modes. The Q factor was not strongly affected, most likely 
because the hole areas were maintained. 

The nanobeam also supports two additional high order modes, shown in 
Figs. 2(c) and (d), at wavelengths 999 nm and 1036 nm, both with Q> Ix 10^, 
which couple to the mechanical mode in Fig. 2(e) with g/2n = 55.2 kHz and 
g/2n = 17.5 kHz, respectively. 

Wavelength conversion 

A cavity optomechanical system supporting two optical resonances at different 
wavelengths and a shared mechanical resonance may be used for performing 
wavelength conversion of classical or quantum optical signals [6, 7, 12, 13]. 
In this photon-phonon translator scheme, a signal tuned to one of the optical 
resonances is transduced to the shared mechanical resonance, and from the lat- 
ter to a second optical resonance at a different wavelength. These conditions 
are achieved in the single nanobeam design above, where the three displayed 
optical resonances couple to the breathing mechanical mode. The wavelength 
separation between input and converted signals is ^ 76 nm, given by the separa- 
tion between the fundamental and third order optical modes. Wider separations 
can in principle be achieved in larger cavities, however, given the requirement 
of small optical and mechanical modal volumes for increased optomechanical 
coupling, g is expected to drop [37]. Ultimately, the achievable spectral sepa- 
ration in single nanobeam designs has an upper bound given by the achievable 
photonic bandgap width, which increases with the refractive index contrast. 
Nanobeams based on Si3N4, with refractive indices n ^ 2.0, are considerably 
more limited than Si nanobeams {n ^ 3.5). 
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Fig. 2. (a) Nanobeam optomechanical crystal, (b)-(d): first, second and 
third order optical resonances, (e) Mechanical resonance. 



4. Slot mode nanobeams 

4.1. Optical cavity design 

We start the design of a Si^N^ optical cavity for operation in the 980 nm band 
by considering the geometry in Fig. 3(a). This is a version of the geometry 
in Fig. 1 in which the mechanical resonator is stripped of its rightmost arm 
and connecting ribs (in what follows, we refer to the remaining arm of the 
mechanical resonator as the 'optical arm'). We first generate Bloch bands for 
photonic crystals with unit cells as highlighted in Fig. 3(a), with optical arm 
widths Wa constant in z. As seen in Fig. 3(b), the fundamental TE Bloch band 
displays bandgaps at the X-point (bands for H^/lVo = 0-25 and H^/Wo = 0.65 
are displayed), and no overlapping bands below the light-line in the spectral 
region near the lowest band-edge. Figure 3(c) shows the squared electric field 
amplitude for the Bloch modes at the Brillouin zone boundary, for the structure 
with Wa = 0.65 • Wq 




Fig. 3. (a) Simplified nanobeam geometry with only the optical arm of 
the mechanical resonator, (b) Photonic bands for TE {E- = at the z = 
plane) mode of the optomechanical crystal unit cell indicated in (a), for 
Wa/Wo = 0.25 and Wa/Wq = 0.65. (c) Normahzed electric field amphtude 
for the fundamental TE mode band at the X point (open circle in (b)). (d) 
Normalized frequency bands for the first four TE photonic Bloch modes 
at the X -point, and (e) normalized effective length for the fundamental TE 
mode, as functions of Wa/Wq ■ 



Figure 3(d) shows the evolution of the first four TE bands at the X-point as 
a function of Wa- The apparent reduction in bandgap (i.e., between the first and 
second bands) width with increasing Wa is due to an effective decrease in the 
filling fraction, or the ratio between air and dielectric regions in the unit cell. 
The larger bandgaps achieved with smaller Wa correspond to a large effective 
reflectivity and thus are better suited for strong spatial confinement. On the 
other hand, a wider Wa is desirable for enhanced modal frequency shifts with 
gap width (i.e., optomechanical coupling). This is evident in Fig. 3(e), where 
Lqm is plotted, assuming |Q ■ n| = 1 in eq. (1), for the fundamental TE photonic 
crystal modes at the F-point, as a function of the Vl^/Wo ratio. A minimized 
optomechanical length can be obtained for Vl^/Wo ~ 0.7, about 30 % lower 
than the maximum value within the plotted range. A trade-off between spatial 
confinement and optomechanical coupling must thus be achieved in a geometry 
where the mechanical beam width Wa is gradually reduced away from the cavity 
center, until a 'nominal' width is reached, for which the reflectivity is highest. 



We thus choose to let Wa vary from 0.65 • Wq at the cavity center to 0.25 ■ Wo at 
the mirrors. 

The photonic bands shown in Fig. 1(a) were obtained with Wy = 0.35 ■ Wq, 
Wx = 0.71 ■ Wo, and t = 0.75 • Wo ( and Wx are the elliptical holes' axes in 
the z and x directions, t is the Si3N4 thickness), parameters that led to a large 
bandgap for Wa = 0.25 ■ Wo- We choose the Si3N4 thickness to be ? = 350 nm 
and the gap to be Wgap = 25 nm, and determine the remaining dimensions by 
regarding the cavity center, where Wo = 0.65 - Wa- A zero-finding routine was 
used to determine the width Wo and the corresponding cavity lattice constant a 
for the lowest X-point TE Bloch mode to be at a wavelength A ~ 980 nm. We 
point out that while such small gap width poses a great fabrication challenge, it 
may in principle be achieved as in [38], where the ?a GPa Si3N4 film stress was 
harnessed to bring two nanobeams together as closely as 40 nm. 

To produce an effective potential well for the fundamental TE mode at the 
X-point, the band edge must shift towards lower frequencies at increasing dis- 
tances from the cavity center. This shift is opposite to what is obtained when 
Wa is modulated in the desired way, and thus we choose to modulate the lattice 
constant along the cavity in order to produce the correct trend. Following the 
procedure described in [28], allowing both beam width Wa and the lattice con- 
stant a to vary quadratically away from the cavity center (Figs. 4(a) and (b)), 
the desired trend for the lowest band edge is achieved (Fig. 4(c)), and an ap- 
proximately linear mirror strength (Fig. 4(d)). The mirror strength corresponds 
to the imaginary part of the Bloch wavenumber for bandgap frequencies at the 
X-point, k= 7r/a{l+iY). From ID first order perturbation theory. 



where £0i^2 are the dielectric and air band edges at the Brillouin zone boundary, 
and cOq is the midgap frequency. Linear mirror strength profiles tend to produce 
optical modes with reduced spatial harmonics above the light line. This leads 
to reduced power leakage into the air, and thus higher quality factors [28]. We 
point out that the quadratic lattice constant and beam width tapers employed 
here, though not optimal, produce sufficiently linear mirror strength profiles for 
high quality factors to be achieved. 

Figure 4(e) shows a cavity mode at a wavelength A = 976 nm, calculated with 
a full vector finite element method, with mechanical beam width and lattice 
constant modulation in Fig. 4(a). The field is strongly concentrated in air gap, 
as in Fig. 3(c), and the quality factor is 2 = 12 x 10^. 





Fig. 4. (a) Photonic lattice constant, (b) optical arm (Wa) over optical beam 
(Wo) width ratio, (c) frequency (at the X-point) of the first two TE photonic 
bands, and (d) photonic mirror strength, as functions of the distance from 
the cavity center (z = 0). (e) Squared electric field amplitude for the pho- 
tonic mode generated with the parameters in (a)-(d). 



4.2. Mechanical resonator design 

The optical resonator design dictates the geometry of the mechanical beam's 
optical arm. The remaining parts of the mechanical resonator (Fig. 1) are de- 
signed around this constraint. While the secondary arm can in principle be cho- 
sen arbitrarily, for simplicity we let it be identical to the optical arm. The lattice 
formed by the connecting ribs are made to follow the optical lattice. This con- 
figuration minimized deterioration of the optical quality factor due to scattering 
at the ribs. Rib widths seemed to have a strong effect on scattering, and had to 
be made small (< 50 nm) to produce sufficiently low Q degradation. 

The mechanical resonator geometry studied here is similar to the optome- 
chanical nanobeam crystal analyzed in [30]. While in this reference a detailed 
account of the phononic crystal mechanical band structure is provided, here 
we focus only on the relevant mechanical bands for enhanced optomechanical 
coupling. We are interested in producing 'breathing' mechanical modes which 
displace the nanobeam arms laterally, allowing the air gap width to change. 

Phononic bands for the crystal shown in Fig. 5(a) are plotted in Fig. 5(b), for 
an arm width Wa = 0.25 ■ Wq- Bands A, B and C were obtained with symmetric 
boundary conditions for the displacement at the / = 0, y = and z' = planes, 
so that no modes with other symmetries are displayed. We refer to ref. [30] 



for a more detailed account of the additional existing bands. The characteristic 
boundary displacements for each band are inset in Fig. 5(b). Near the F-point, 
the mechanical modes of band C display the desired lateral displacement pat- 
tern, necessary for the formation of breathing resonances. 

Modes from bands A and B stem from flexural modes of the connecting 
ribs - evident in the inset displacement plots- and display negligible lateral dis- 
placement. Band C has a minimum at kj^ = 0, which creates the conditions for 
the formation of a breathing mode phononic resonance. Indeed, a phononic 
cavity naturally arises from the quadratic spatial arm width (Wa) variation, as 
evidenced in Fig. 5(c), where the F-point edge of band C is plotted as a func- 
tion of Wa/W^o- As the arm width decreases from the cavity center towards the 
mirror regions, the band edge moves towards higher frequencies, placing the 
resonance frequency within the breathing mode phononic bandgap. Localized 
breathing modes as plotted in Figs. 6(a) and (b) result, with frequencies that 
decrease with increasing nanobeam widths Wm, as shown in Fig. 6(c). It is also 
apparent that displacement profiles can also change considerably with Wm- This 
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Fig. 5. (a) Phononic crystal geometry, (b) Mechanical Bloch mode band 
structure for a phononic crystal with Wa /Wq = 0.25. The inset shapes show 
a top view of the unit cell, with exaggerated boundary displacements for 
modes on bands A, B and C. (c) Evolution of the F-point eigenfrequency 
for band C, as a function of the ratio between the optical arm (Wa) and 
optical beam (Wq) widths. 



is quantified in Fig. 6(d), where the normalized displacement 

Z,,.Ji^, (3) 

J dA\Q-n\ 

is plotted as a function of Wm/Wq- The integral is performed over the optical 
arm's surface facing the optical beam, and so is a measure of the uniformity 
of the optical arm's displacement towards the optical beam. The evolution of 
the modal displacement pattern also results in a variation of the motional mass, 
as shown in Fig. 6(e) 

Optomechanical coupling 

The optomechanical coupling goM and coupling rate g between the optical res- 
onance calculated in Section 4. 1 and the breathing mode mechanical resonances 
discussed in Section 4.2 are plotted in Fig. 6(f), as a function of Wa/Wm- The 
maximum ^qm and g values achieved for the Wm/Wq range shown are > 2x 
the value obtained with the optimized single nanobeam of Section 3, despite 
the larger motional masses (Fig. 6(e)). The strong influence of the gap width 
on the optomechanical coupling is evident in the fact that both g and goM ap- 
proximately follow the trend of D^. In fact, the effect of increasing gap widths 
on the resonance wavelength is plotted in Fig. 6(g), displaying a fast blue shift 
for decreasing gaps. Because the cavity was initially designed to have a 25 nm 
gap, the quality factor decreases as this parameter is changed, and so do the 
optomechanical coupling goM and the coupling rate g, despite the increase in 
the optical frequency {g oc x^pf oc /„, x^pf = \Jh/2m(0m is the zero point fluc- 
tuation amplitude for mechanical displacement). For any fixed gap width, how- 
ever, the hole lattice can be adjusted to yield a high quality factor. For instance, 
changing the the lattice constant profile in Fig. 4(a) so it varies (quadratically) 
from 325 nm at x = to 365 nm at the cavity edge gives an optical mode at 
wavelength A = 986.4 nm with 2 = 2.4 x 10^ and goM = 180 GHz/nm. In con- 
trast, for the 25 nm gap design of Fig. 6, a 40 nm gap yields X = 965.6 nm, 
2 < 5 X 10^ and gOM = 150 GHz/nm. This serves to show that our double 
nanobeam design affords great flexibility towards optimized performance un- 
der constraints such as resonance wavelength, quality factor or gap width. 

4.3. Effect of index contrast 

Equation 1 suggests that the higher index contrasts generally yield higher op- 
tomechanical coupling. In addition, the electric field concentration in the slot 
region (oc (An)^ [14]) is expected to add a significant contribution to goM 
and g. Indeed, a design for operation at the 1550 nm band based on sili- 
con (nsi ~ 3.48), produced with the same procedure as above, yielded an op- 
tomechanical coupling rate gjln ^ 900 kHz for a mechanical resonance at 
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Fig. 6. (a)Fundamental optical gap resonance, (b) and (c): Mechanical 
beam breathing modes for(b) Wm/Wq = 2.6 and (c) Wm/Wq =2.8. (d) 
Mechanical frequency, (e) normalized displacement D^, (f) motional mass 
m and (g) g and ^om for the fundamental breathing mode as a function 
of Wm/Wq- As a function of gap width: (h) fundamental optical resonance 
wavelength, (i) quality factor and (j) optomechanical couping (^om) and 
coupling rates (g) with the fundamental mechanical breathing mode for 
Wm/Wq = 2.8 of (b). Note that the geometry that produced (h)-(j) was 
optimized for a 25 nm gap. Improved values of gou, g and Q than plotted 
can in be obtained with proper optimization for different gap widths (see 
text). 



fm = 1.38 GHz (an optimized single nanobeam design with similar geomet- 
rical parameters displayed g/27Z ^ 400 kHz for a mechanical resonance at 
fm = 4.3 GHz). Table 1 shows the relevant optomechanical quantities for this 
design. Here, the ID lattice constant was allowed to vary quadratically from 
360 nm to 385 nm, and the optical beam width varied between from 325 nm at 
the cavity center to 150 nm. 



Table 1. Silicon Based Optical Resonator Optomechanical Crystal Param- 



eters. A slot width of 10 nm was assumed. 



K (nm) 


/«,(GHz) 


<2x 10^ 


gOM (GHz/nm) 


gjlTl (kHz) 


LoM ( Aim) 


m (pg) 


1540.4 


1.38 


1.5 


485 


882 


0.4 


1.8 



5. Applications 

5.1. Wide spectral separation wavelength conversion 

While a 'standard' nanobeam geometry such as that studied in Section 3 al- 
ready offers the necessary conditions for wavelength conversion, the spectral 
separation between optical resonances of different orders is limited to, at best, 
70 nm. Even in silicon nanobeam designs, where the higher index contrast 
leads to a considerably wider spacing between resonances, a maximum sep- 
aration of ^ 100 nm is achievable [12]. As discussed in Section 3, the achiev- 
able wavelength separation in the single nanobeam design is ultimately given 
by the nanobeam photonic bandgap, which can be limited even for high con- 
trasts such as in the Si case. The ability to perform wavelength conversion over 
considerably wider wavelength separations may be desirable for applications 
in which classical or quantum information transmitted over optical fibers at 
telecom wavelengths is originally produced at completely separate wavelength 
bands, for instance in the near-infrared or visible range. Such an achievement 
would be enabling for proposed hybrid quantum optical networks, whose nodes 
may be widely different from each other (trapped atoms [39], quantum dots [40], 
organic molecules [41], etc). 

The physical separation between optical and mechanical resonances afforded 
through our double nanobeam optomechanical approach allows the formation 
of optomechanical resonators with widely spaced, spatially separate and semi- 
independent optical resonances with reasonable optomechanical coupling to 
the same mechanical resonance. Such a resonator, schematically depicted in 
Fig. 7(a), consists of a central mechanical nanobeam laterally sandwiched by 
two optical nanobeams. Air gaps on the two sides of the mechanical resonator 
allow a single breathing mode resonance to couple efficiently to optical nanoslot 
modes on both sides. We demonstrate the feasibility of this scheme via an ex- 
ample in which optical resonances at wavelengths of 980 nm and 1310 nm are 
supported. Figures 7(b) and (d) respectively show the 1310 nm and 980 nm res- 
onances of a double nanobeam resonator that share the mechanical resonance 
at /,„ = 1.38 GHz in Fig. 7(c). Both optical resonances display quality fac- 
tors in excess of 10^, which allows operation in the resolved sideband regime. 
The 980 nm wavelength optical cavity has exactly the same geometry param- 
eters as that in Section 4.1. The 1310 nm cavity was designed with the same 



procedure as described in Section 4.1. Achieving high optical quality factors 
for both cavities required significant modifications to the original mechanical 
nanobeam geometry presented in Section 4.2. First, the nanobeam width Wm 
had to be increased to 4 ■ Wq, to minimize the interaction between the two op- 
tical resonators. This resulted in a lower resonance frequency, /„, = 1.38 GHz. 
In addition, the 1310 nm resonator's optical arm was made constant with 
= 0.25 ■ Wo, to prevent coupling of the confined 980 nm resonance to leaky 
resonances of the 1310 nm beam. This caused the optomechanical coupling 
gOM for the 1310 nm resonator to be ~ 40 % lower than that obtained with 
modulation of the optical arm (goM ~ 1 12 GHz/nm). It is worth noting on the 
other hand that, despite the uniform arm width, goM can still be reasonably 
high. This parameter, together with the optomechanical coupling rate g, effec- 
tive optomechanical length Lqm and motional mass m for the two resonators 
are given in Table 2. 



Table 2. Double optical Resonator Optomechanical Crystal Parameters. 



Both s] 


ots are assumed to be 25 nm wide. 


Xo (nm) 


Qx 10^ 


gOM (GHz/nm) 


g/2;r (kHz) 


Lqm ( jUm) 


m(pg) 


978 


4.4 


179 


322 


1.7 


1.92 


1318 


4.6 


67 


122 


3.4 


1.85 



As discussed in the Appendix, the conversion efficiency achievable in the 
scheme proposed in ref. [6] and demonstrated in ref. [12] depends on the co- 
operativities Ci^2 °^ gf 2^c;i,2 of optical cavities 1 and 2 (nc;i^2 is the average 
photon population in each cavity), and unity conversion requires Ci = C2 ^ 1. 
The coupling rates calculated above are approximately a factor of three times 
smaller than those used in ref. [12]. These smaller coupling rates can be com- 
pensated by increasing the number of photons in the cavity (by a factor of nine, 
assuming similar optical and mechanical decay rates are achievable), which 
should be feasible in the Si3N4 system. We also note that the disparate coupling 
rates for the two cavities is not expected to influence the ability to achieve high 
efficiencies, as such differences may be compensated by proper balancing of 
the intra-cavity photon populations ric- 1,2- 

5.2. Strong optomechanical coupling 

The observation of quantum behavior in cavity optomechanics relies on the 
possibility of strong coupling between optical and mechanical resonances. Ulti- 
mately, a regime is sought in which interactions at the single photon and phonon 
levels are observable. This single-photon strong coupling regime, in which a 
single phonon is able to shift an optical resonance by an extent comparable to 
the latter's linewidth, is characterized by a coupling rate g comparable to both 
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Fig. 7. (a) Double optical cavity optomechanical crystal geometry, (b) 
1310 nm optical mode, (c) Mechanical resonator breathing mode, (c) 
980 nm optical mode. Green arrows indicate simultaneous coupling be- 
tween the optical resonances and the mechanical mode. 



I 



the optical resonance decay rate k and the mechanical frequency Q.. Achiev- 
ing single photon strong coupling at optical frequencies is challenging due to 
high optical losses. In [10], a scheme was proposed for producing a strongly en- 
hanced, effective optomechanical coupling for quantum non-demolition (QND) 
photon detection and phonon number readout measurements. Such a QND 
measurement would extract information from the quantum system without dis- 
turbing its quantum state. In the setup of [10], an optomechanical crystal geom- 
etry supports two optical resonances split by a frequency 27, comparable to the 
mechanical frequency Q. of an interacting mechanical resonance. The coherent 
interaction between photons and phonons is characterized by an effective cou- 
pling rate ^q/ 5Q., where 5Q. = 2J — Q. and is the bare optomechanical cou- 
pling rate. In the limit 27 — Q, the optical frequency shift is g^/ 6Q.nh, where 
nti is the phonon number,which may be sufficiently large to produce a phonon 
number readout. 
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Fig. 8. Optical mode splitting for the symmetric double optical beam res- 
onator in the inset. The dotted line in the inset indicates the symmetry plane 
that defines splits the original optical resonance into symmetric (blue) and 
anti-symmetric (red) modes. Parameters for the optical and mechanical res- 
onators are as in Section 4, with Wm/Wq = 2.8 



The flexibility of our double nanobeam geometry allows us to accommodate 
the optical mode splittings 2J ^ Q. required in the scheme proposed in [10]. 
As shown in the inset of Fig. 8, an optomechanical resonator is formed by two 
identical optical beams sandwiching a mechanical resonator. The two individual 
optical resonances couple evanescently, forming a symmetric / anti-symmetric 
doublet whose spectral separation can be controlled by the mechanical beam 
width Wm- The modal splitting 27, shown in Fig. 8 for the same cavity param- 
eters as in Section 3, indeed decreases with increasing Wm, as the spatial over- 
lap between the individual optical resonances diminishes. The main breathing 
mode mechanical frequency also decreases with Wm, however remaining within 
the GHz range even for the largest mechanical beam widths plotted, so that 
the regime 5Q. = 27 — ^2 — )■ can be achieved. The coupling between the two 
modes can also be tuned by creating optical cavities with coinciding individual 
resonances, but different gap widths. 

5.3. Microwave to optical photon conversion 

State-of-the art superconducting microwave circuits can be designed to provide 
controllable, coherent interactions between superconducting qubits and high 
quality microwave cavity resonators, and provide favorable conditions for the 
creation and manipulation of photonic states at microwave frequencies. In the 
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Fig. 9. Schematic for possible optomechanical microwave to optical wave- 
length converter. The RF cavity resonance frequency is modulated by the 
gap capacitance Cg, which depends on the gap width, and therefore on the 
displacement of the mechanical resonator. The latter is coupled to the slot 
optical mode formed with the optical beam. 

context of a quantum network [39] where the nodes are connected via photonic 
links, however, low loss transmission of microwave photons over long distances 
becomes a challenge, as superconducting transmission lines would be required. 
The ability to convert between microwave and optical photons which can be 
transmitted via optical fibers thus constitutes an appealing solution to this prob- 
lem. 

The physical separation between optical and mechanical resonances in our 
geometry creates opportunities for the generation of microwave-to-optical fre- 
quency signal transducers. One possible geometry is illustrated in Fig. 9, where 
the far arm of the mechanical resonator incorporates a metallic strip, and, to- 
gether with a neighboring ground plane, forms a capacitor that integrates a su- 
perconducting, planar circuit microwave resonator, similar to that demonstrated 
in [42, 43]. The displacement of the mechanical beam affects the cavity reso- 
nance by causing a variation in the circuit capacitance. Coupling between the 
mechanical displacement and cavity capacitance is stronger for smaller gaps 
between the metallized beam and the ground plane. At the same time, the dis- 
tance between the metallized mechanical beam arm and the optical resonator 
must be sufficient to prevent significant deterioration of the quality factor of 
the optical cavity, while keeping the superconducting microwave resonator suf- 
ficiently separated from strong optical fields is necessary to avoid heating and 



degradation of its performance. As shown in Section 3, this separation can be 
accomplished through increased mechanical beam widths. While this leads to 
decreasing mechanical frequencies for the fundamental breathing mode, as seen 
in the previous sections, mechanical frequencies in the GHz range are achiev- 
able. 

As is the case with optical-to-optical conversion, microwave-to-optical con- 
version under the protocol described in ref . [6] requires C^wave = Copticai ^ 1 , 
that is, the cooperativity of the microwave and optical resonators (see Ap- 
pendix) should be large and matched to each other. We anticipate that the design 
of the optical resonator will be similar to that described previously, while a full 
design of the microwave cavity is beyond the scope of this paper. Qualitatively, 
we note the conceptual similarity between the geometry shown in Fig. 9 and 
that studied experimentally in ref. [43], where a suspended Al nanobeam sepa- 
rated from a ground plane by a gap ^ 20 nm was employed to perform vibration 
detection near the quantum limit, and for which subsequent experiments have 
shown radiation-pressure driven phenomena like electromagnetically-induced 
transparency [44]. The coupling between the microwave cavity and mechanical 
resonator is ^^wave ~ {dCg/dwgap), where Cg is the coupling capacitance and 
Wgap is the gap between the two elements. Reaching the regime C^wave 2> 1 with 
a microwave intracavity photon number consistent with recent experiments (up 
to 10^ photons were used in ref. [44]) will thus require small gaps (10 nm gaps 
were reported in ref. [43]), along with low dissipation for both the microwave 
cavity and the mechanical resonator. Investigation of such suitable microwave 
cavity geometries is in progress. 

6. Conclusion 

We have presented a design methodology and analysis of an optomechanical 
crystal comprised of two laterally coupled ID nanobeams which support spa- 
tially separate optical and mechanical resonances. This spatial separation en- 
ables independent design and flexible optimization of the optics and mechanics 
of the system, and can be of particular importance in applications requiring the 
coupling of multiple electromagnetic modes to a single localized mechanical 
resonance. In addition, the small gap (^25 nm) separating the two nanobeams 
gives rise to a slot optical mode effect that enables a large zero-point optome- 
chanical coupling strength to be achieved, even for nanobeam materials with 
small refractive index contrasts such as Si3N4. Our design predicts optical qual- 
ity factors above 10^ and mechanical frequencies above 1 GHz for nanobeams 
in both the Si3N4 and Si systems. While the optomechanical coupling strengths 
predicted here are moderately superior to those calculated and experimentally 
determined from single nanobeam geometries, optical and mechanical reso- 
nances can be designed with considerably more independence and control. 



The large predicted optomechanical coupling strengths to GHz mechanical 
oscillators in Si3N4 are particularly attractive, as this material displays a broad 
optical transparency window, which allows operation throughout the visible and 
near-infrared. The material also offers a low intrinsic mechanical dissipation 
rate, and does not exhibit the two-photon absorption and subsequent free-carrier 
absorption and dispersion observed in silicon [27], potentially allowing higher 
powers to be employed and increasing the range of achievable pump-enhanced 
optomechanical coupling values. All of these features provide good prospects 
for the realization of radiation-pressure mediated photon-phonon translation 
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Appendix - Optomechanical wavelength conversion 

The wavelength conversion scheme proposed in Sections 3 and 5.1 involve two 
optical resonances at frequencies cOo.i and C0o,2, both of which are coupled to a 
mechanical resonance at frequency (y,„. Both optical cavities are prepared with 
strong pumps, red-detuned from the resonance centers by the mechanical mode 
frequency. For each optical resonance, this preparation gives rise to an effective 
optomechanical interaction/or signals at the cavity center described by a beam- 
splitter type interaction Hamiltonian //int,eff = hG {a^b-\-db^^, where a and b 
are the destruction operators for cavity photons at the center frequency and 
phonons in the mechanical resonator. The parameter G is the pump-enhanced 
optomechanical coupling, G — g\(X,^^\, where is the square root of the steady- 
state photon population at the pump frequency, and g the bare-cavity optome- 
chanical coupling rate. It is apparent from the effective Hamiltonian that the 
interaction allows the exchange of signal photons and phonons with a rate dic- 
tated by the intensity of the pump beam. In the wavelength conversion process, 
a signal photon at frequency (Oi is injected into the (prepared) cavity 1, and con- 
verted to a cavity phonon. The process is reversible, such that the newly created 
phonon may be subsequently transduced to a photon in (also prepared) cavity 
2, at frequency (Oo2- This process can be described classically via the matrix 
equation [12] 



[6]. 
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(A-1) 



obtained from the Heisenberg equations for the bosonic operators ai 2 (for 
photons in cavities 1 and 2) and b. To arrive at eq. (A-1), we replaced 
«i,2 — ^ ai^2^^'^'"'-''^'^^'^\ j8 — )■ /3e^'^^ where A is the detuning of the optical 
signal from the cavity center frequency C!)c;i,2> and the rotating wave approx- 
imation was used. Noise sources were completely disregarded. In eq (A-1), 
K"i,2 = 'Qx;i,2 + K';;i,2 are the decay rates for optical modes 1 and 2, comprised 
of intrinsic (k,, associated with e.g., radiative losses), and extrinsic (Kgx, associ- 
ated with coupling to external access channels) components; jj is the intrinsic 
decay rate for the mechanical mode; Win and is the input optical field. Finally, it 
was assumed that the detuning between the pump and signal beams (at the in- 
put and output optical cavities) was equal to the mechanical frequency, A = cOm- 



The field output from cavity 2 is simply aout;2 = v K'ex,2/2a2, so eq.(A-l) can 
be solved to give 

«out;2 = VnVni- «in,l , (A-2) 

7/2 

where 7 = 7om;1 + roM;2 + Yi and 7om;1,2 = AGI ^/ ^1,2 (roM;i,2 are the opti- 
cal spring contributions to the total mechanical mode damping 7), and 771,2 = 
'Qx;1,2/2k"i 2.The conversion efficiency, then, is 



0Jout,2 



0^in,l 



4C1C2 



where C\2 = 4Gj 2/ K'i,27 are the cooperativities for optical cavities 1 and 2. It 
is apparent that, for T] — )■ 1, Ci = C2 ^ 1 and r/i 2 — ?• 1 are sufficient. A more 
complete derivation of the wavelength process,including treatment of noise, is 
given in [12]. 



